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Abstract 

Considering the high-energy hmit of the QCD gluon distribution inside a 
nucleus, we calculate the proton-nucleus total inelastic cross section using 
a simplified dipole model. We show that, if gluon saturation occurs in the 
nuclear surface region, the total cross section of proton-nucleus collisions 
increases more rapidly as a function of the incident energy compared to that 
of a Glauber-type estimate. We discuss the implications of this with respect 
to recent ultra-high-energy cosmic ray experiments. 
Keywords: 
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1. Introduction 

In the ultra-high-energy domain, the mechanism of inelastic hadronic 
collisions is dominated by the contribution from small-x gluons. This is the 
basic reason why the total or inelastic proton-proton collision cross section 



Email addresses: licinioOif .uf rj .br (L. Portugal), tkodamiaOif . uf r j .br (T. 

Kodama) 

Preprint submitted to Nuclear Physics A March 11, 2010 



increases as a function of the incident energy. A simple way to see this is 
to use the eikonal expression for the reaction cross section in the impact 
parameter representation as 

ar = J d2b[l-exp(-2x(b,s))], (1) 

where the eikonal x(t>, s) counts essentially the total number of possible 
scattering centers of the constituents inside the target "seen" by the projectile 
passing through the target along a straight line with the impact parameter 
b, and with center-of-mass energy ^/s. 

If the target is thick enough, the eikonal function is much larger than 
unity ( x(b, s) ^ 1) in the central region, but it falls down to zero in the 
surface region. The integrand of Eq. ([T]) keeps a value almost equal to unity 
while x(b, s) ^ 1, and falls down quickly to zero near the surface. From 
this we may determine an effective radius bi/2 such that the integral ([1]) is 
approximately given by 

ar ~ vr6i/2- (2) 

Here bi/2 might be estimated by fixing the value of x0^i/2,s), for example, 

as 

X(bi/2,s) = ln2, (3) 

so that the effective radius depends on the energy bi/2 = bi/2 (y/s). Let us 
assume that the eikonal function factorizes in the form 

x{h,s) = P{h)N{V~s), (4) 

where -P(b) is the probability distribution of the scattering center (e.g., par- 
tons) in the transverse plane and N{x) is the number of partons which can 
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interact for a given ^/s. Suppose that -P(b) is a two-dimensional Gaussian 
distribution of width i?, 

and that, for large -y/i, the number of partons increase with ^/s as 

A^(a;) oc v^". (6) 

Then, we have 

52^2 (v^) = aR^ In + Const. (7) 

In such a situation, the reaction cross section increases as function of the 
incident energy for very large ^/s as 

ar ^ avri?^ In \/s. (8) 

If the edge of the distribution has an exponential tail instead of a Gaus- 
sian distribution, a similar argument will show that the cross section would 
increase as 

CFr ~ Const X (in a/s)^ . (9) 

The important point of this simple argument is that the rate of increase 
is related to the diffuseness R of the probability distribution of the scat- 
tering center. The more diffuse the surface thickness is, the more quickly 
the reaction cross section increases as function of the incident energy. The 
possibility of such a mechanism, not only in the proton-proton but also in 
nucleus-nucleus cross section, was suggested many years ago |2|. 

In this paper we explore the idea of [l|,l2| in the language of the QCD gluon 
saturation mechanism for the proton- nucleus reaction. We show that, if the 



gluon distribution becomes saturated at some energy scale inside the nuclear 
surface region, then the reaction cross section of proton-nucleus collisions 
starts to increase very quickly and eventually overcomes the values estimated 
by the usual Glauber type of calculation. 

Applying a simple effective dipole model for the reaction mechanism, 
we find that such an energy scale is of the order of 10^'' — lO^^eK Above 
this energy scale, the behavior of the proton-nucleus cross section begins to 
change. We suggest that such an energy dependence of the proton-nucleus 
cross section may be observed in terms of the quantity called (Xmax) in the 
air showers of ultra-high-energy cosmic rays. Using a very simple toy model 
estimate of (Xmax), we show that our calculated values of the proton- nucleus 
reaction cross section are consistent withjthe recently observed (X^ax) by the 
Pierre Auger Observatory experiments pi] for incident protons at ultra-high 
energies. 



2. Effective Dipole Model for the Proton-Proton Cross Section 



To explore the above idea, we need the total cross section of proton-proton 
as function of incident energy which permit to extrapolate to the ultra-high 
energy region. For this purpose, we introduce an extremely simplified version 
of the dipole saturation model for the purpose of illustratirig our idea here. 
The original dipole model was first introduced by Mueller ^ and extended 
to the impact parameter representation by Kowalski and Teaney jsj. Re- 
cently this model was applied to fit diffractive structure functions in electron- 



71, Isl]. An application to proton-proton processes can be 



Nucleus collisions 

n 

found in ^], where the incident proton is represented by dipoles associated 
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with the valence quarks, having a probabihty distribution for different sizes 
specified by the proton wave function. 

In this work, for simphcity, we assume that the proton is described by an 
effective dipole of a given size Re, where we fix the dipole size -Rd by the 
average dipole radius 

Rl = j d?rTrl\^p{rT)\\ (10) 
where \l/p(rT) is the proton dipole wave function 10| in the transverse plane. 



In this case, the gluon saturation model leads to a reaction cross section, 

/•OO 

Or (v^) = 27r / hdh (1 - e-2>^) , (11) 
Jo 

where the eikonal for the proton-proton reaction can be written as 
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X(b, s) = —Ri,asiQ')xgix, Q^)Tp(b). (12) 



' c 
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In this expression, the quantity xg{x, Q ) is the parton distribution function 
(PDF), and represents the gluon density x— distribution in the target at a 
scale Q^, and asiQ"^) is the strong coupling defined at this scale. Following 



the steps of 



8| we use 



Q" = ^ + Ql (13) 



D 



5 



with Nf = 3 and Aqcd = 0.192 GeV. 

For the profile function Tp(b) of a proton, we consider the two different 
cases described in below. 



2.1. Gaussian profile function 

Suppose that the transverse probabihty distribution of gluons (profile 
function) inside the proton is of a Gaussian type: 

T,(b) = ^axp(=|). (16) 

Writing Eq. flTT]) as 

ae 



(Jr (v^) = 27ri?p j xdx {l — e~ 
Jo " 



where 



= ttRH^e -Chi (a) + Shl{a) + \n{a)} 

= 7iRlF{a), (17) 



^'''as{Q')xg{x,Q'), (18) 



jE is the Euler constant, Chi (a) and Shl{a) are hyperbolic-cosine and 
hyperbolic-sine integral functions. The function 

F{a) = {'jE - Chi (a) + Shl{a) + In (a)} (19) 

is a monotonic increasing function of a. From these expressions, we have 
a oc a/s" and F (a) ~ In (a) for very large i/i, so that the reaction cross 
section behaves at large incident energies as 

ar ~ a ttRI In ( v^s) + Const. (20) 

as expected from the simple argument presented in the Introduction. 



2.2. Profile function with exponential tail 

Another important possibility is that the spatial gluon distribution lias an 
exponential tail instead of a Gaussian, such as in a Woods-Saxon distribution. 
In this case, the asymptotic increase of the cross section as function of energy 
is expected as log^ (v^) and not log {^/s) , as discussed in the Introduction. 
However, since the Woods-Saxon distribution has 2 parameters (half-density 
radius R1/2 and surface thickness d), wc use the following profile function 
which has exponential tail, for simplicity. 



where Ct is the Catalan number, 

{2i + l)\ 



°° (-lY 

'^'• = E74tW-»-91«. (22) 



i=0 

The width parameter Rp is related to the mean-square of the impact param- 
eter (6^) as 

{b') = i^Rl (23) 

and consequently to the radius parameter of the Gaussian case as Rp — 
y/vRp, where 

fn°° I coshx dx , , 

^ = f °o . , - 6-478. 24 
Jq xj coshx dx 

Then: 



xdx h-e""^^ 



cosh{x) 



= 7rRlx2uI{a), (25) 
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where we defined the integral 

xdx (^1 -e"5%/'°'^(")j , (26) 

where a is given in Eq.( IT8ll . 

We can examine the asymptotic form of this integral for large a as follows. 
Changing variables, t = a/ cosh a; with a = a/AuCr, the integral above can 
be rewritten as 

We then separate the integral into two parts: 

where /3 is some finite constant but sufficiently larger than unity so that 
<^ 1. For large a ^ /3, the dominant asymptotic contribution comes from 
the second integral. In this integral, the second term e~* may be neglected 
when compared to unity, and we have 

Therefore, the total reaction cross section at large ^/s behaves as 

ar ~ 2z/7ri?J In^ [^s) + O (in v^) , (30) 

as discussed above. 

In Fig.l, we show the behavior of the functions F (a) and 2vl (a). Notice 
that for a < 0.5 the two functions almost coincide, but for a > 2 they behave 
as ~ In (a) and ~ In [of , respectively. 
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2.3. Fits to the proton-proton cross section 

Naturally, our simple effective dipole description will not work well for 
the low-energy region. However, the objective of the present work is to show 
the effect of possible gluon saturation inside the nuclear surface region at 
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high energies, we just readjust slightly the parameters determined in 
fit t 



to fit the er ierg y dependence of the proton-proton reaction cross section 
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14 



, ll5|] only for y/s > 100 GeV . Note that the cross sections for 
^/s ~ AOTeV are from the cosmic ray data extracted from the proton-light 
nuclei interactions. We can obtain reasonable fits using both Gaussian and 
hyperbolic-secant profile functions at higher energies, as seen in Fig.2. Circles 
are experimental data, crosses represent the result of the Gaussian profile 
function, and open squares are the result of the 1/ cosh profile function. In 
this energy region, both curves are similar but, as expected, the cross section 
for the 1 / cosh profile function shows a more rapid increase in energy at higher 
energies than the Gaussian profile case. For these calculations we have used 
the PDF xg{x,Q'^) from the GRV98 collaboration |l6|. 



3. Proton-Nucleus cross section 

3.1. Independent Nucleon Glauber Picture (INGP) 

For a high energy proton-nucleus collision, we may calculate the reaction 
cross section ap^A as a superposition of independent nucleon-nucleon colli- 
sions in the Glauber approach. Hereafter, this picture is referred to as INGP. 
In this picture, we have the well-known formula 

ap+A = j d'^h (l - e-^^-(^)'^-(^)) , (31) 
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where Upp is the total cross section of the proton-proton coUision, and 
Tat (b) = dz PA(b,-2) is the transverse probabihty distribution function 
of a nucleon inside the target nucleus. 

The derivation of this formula can be found in many places. Here, for later 
comparison, we summarize it. First, assume that the partons are confined in 
each nucleon of the target nucleus and, for each coUisional event, nucleons are 
independently distributed inside the target nucleus. In this vision, for a one 
collisional event the incident proton "sees" the target nucleus as a collection 
of nucleons whose center-of-mass positions are specified as 



in the transverse plane. For this event, we have the eikonal expression of the 
total reaction cross section 



where Xvp (^5 V^) eikonal for proton-proton collisions of impact pa- 

rameter b. 

To calculate the proton-nucleus reaction cross section, we should take the 
average over all events. Assuming that the target nucleus is an ensemble 
of independent nucleons whose single particle distribution is Pa {t) , we can 
write the average over events as 



{bi, b2, .., b^} 



(32) 




(33) 




(34) 



where 
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When -Pa(?) is a slowly varying function compared to the nucleon size, we 
can write 




We identify the term 

J d%' (l - e-^-(''-^)) = (yi) , (37) 

so that 

F(b, V^) (v^i) T, (b) , (38) 

where 

Tn (h) = I dz P{h,z) (39) 

is the transverse distribution function of a nucleon inside the nucleus. 

Eq.( l38|) becomes negative for appT^ > 1 while in Eq. (!35|) it is defined as 
positive definite. Considering the shadowing effect, we should replace Eq. (!38|) 
by an eikonalized expression, 

^ _i. g-T'.(b)<7j,j,(Vi)_ ^^Q^ 

From this, we have 

(ap+A) = J d'h (l - e-^^-(b)-pp(Vi) j ^ (41) 
which is nothing but Eq.fj3T|). 
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Since 

(v^) - In (v^) , (or In^ (7^)), (42) 

we conclude that the INGP leads to a very weak energy dependence of the 
cross section 

~lnln(v^) (43) 
for large y/s, and Eq. fllTl) gives essentially the geometric cross section of order 

We must remind that, in Eq. pT]) the effects of finite size of nucleon and 
of the fluctuation were not included. 

3.2. Gluon saturation in the nuclear surface (GSNS) 

In contrast to the approach above, we may consider the proton-nucleus 
collision process in terms of the gluon saturation inside the whole nucleus. 
When we go to sufficiently large energies, gluons of bounded nucleons inside 
a nucleus should start to superimpose and eventually fill up the nucleus as 
a whole. In this regime, we should then use the dipole model with gluon 
saturation inside the nucleus to calculate the total cross section for proton- 
nucleus collision. Hereafter, such a scenario is referred to as GSNS. 

The proton nucleus cross section is then 

(TpA [y/s) = 27r / hdh (1 - e-2x) , (44) 

where the eikonal for the proton-nucleus reaction can be written as 

X(b, s) = ^Rla,iQ^)xgix, Q')Tr,{h). (45) 

In this regime, as we discussed before the total reaction cross section in- 
creases as In {^/s) or In (s)^ , depending on the form of gluon profile function 
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near the nuclear surface which is much quicker than In In {^/s) of the indepen- 
dent nucleon picture. Therefore, the total cross section for proton-nucleus is 
eventually dominated by the gluon saturation process inside the nucleus. 

It is interesting to investigate in what energy scale such phenomena will 
happen and what physical parameters determine such a scale. In the inde- 
pendent nucleon picture, the energy dependence of the total cross section is 
very weak and stays more or less at the order of the geometrical cross section, 
nRj^. Therefore, we conclude that the crossover energy scale should happen 
when F (a) ^ 1 or 2vl {a) ~ 1 from Eqs. fll7ll25p . From Fig.l, we find that 
this occurs at 

a ~ 1, (46) 
where a is given in Eq.f fTSj) . For large ^/s, we may approximate the PDF as 

xg{^s,Q)^^{Q) (47) 

The above condition gives an estimate of the energy scale for the crossover: 

(48) 

From this expression, we conclude that the energy scale for the gluon satura- 
tion scenario inside the target nucleus is determined essentially by the ratio 
Rn/Rd- The smaller the ratio is, the smaller the energy scale becomes. 

4. Ultra high energy cosmic rays 

To see in practice what is the crossover energy scale for the gluon satura- 
tion inside a nucleus, we compare the calculated total reaction cross sections 
using the INGP and the GNSN visions. We take typical air nuclei of average 
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Crossover 



D_ 



c^s{Q')aQ) 



{A) = 14.5 where = 1.1 A^^^ fm. We also calculate the proton-air nuclei 
cross sections for the following 3 cases of different nuclear profile functions: 



• Gaussian, Eq.f lT6|) substituting Rp by Rn- 



• Hyperbolic secant, Eq. fl2T]) substituting Rp by Rn 



• z-Integrated Woods-Saxon 




oo 



1 



dz 



1 + exp { {y/WT^ - Rn) /a 



} 



(49) 



— oo 



with a = 0.5 fm and normalization factor 



) 



2 



) 



(50) 



In Fig. 3 we compare the energy dependences of proton-air collision cross 
sections calculated for various situations: One in the INGP (■) , and other 3 
cases of the GNSN with Gaussian nuclear profile (+), GNSN with hyperbolic 
secant profile (x) and the 2— integrated Woods-Saxon profile (□). The INGP 
gives almost the same cross section for all the profile function so that only 
one line is shown. 

In these calculations, the parameters of the model {Rn, Rp, C, Ql ) are 
those fitted to the proton-proton cross section using the Gaussian profile 
(see Table I). In order to calculate the cross section in the ultra-high-energy 
cosmic ray domain (> lO^^el^) , we have to extrapolate the PDF for small x 
values. For the Gaussian profile we have = 18.8 GeV^ and, consequently, 
we fit the PDF as 



xg (x) = 6.56x 



-0.271 



(51) 
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For the 1/ cosh case we have = 17.2 GeV^ and the corresponding PDF 
becomes 



Naturally, these extrapolations contain ambiguities. However, as we are in- 
terested in obtaining a schematic description of the proton-proton cross sec- 
tion, these ambiguities can be absorbed in the final fitting procedure of pa- 
rameters of the effective dipole model so that they will not affect our general 
conclusion. 

As expected, the GNSN scenario gives a rapid increase of the proton- 
nucleus cross section as a function of the incident energy and eventually 
overcomes the value of the INGP. It is interesting to note that all of the cross 
sections for different profile functions of GNSN cross the INGP estimates at 
the energy scale of 10^^ — lO^^eV. This is because, once the proton cross 
section is fitted, the values of a defined in Eq.f jTSj) are not much different. 

As mentioned before, the use of the 1 / cosh profile function for the gluon 
distribution in a proton gives a more rapid increase of the p-p cross section at 
larger energies compared to the Gaussian profile function. When we use this 
type of fit to the proton-proton cross section, it naturally leads to a larger 
proton-nucleus cross section for the INGP, as shown in Fig. 4. However, the 
proton- nucleus cross section calculated in the GSNS scenario stays invariant. 
This is because, once the PDF and effective dipole parameter are determined, 
the proton-nucleus cross section in this scenario depends only on the gluon 
distribution inside the target nucleus. 

One direct consequence of such effects should reflect in the behavior of 
the observable (X^ax) , essentially the normalized depth of the position of 




(52) 
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maximum luminosity of an air-shower in the atmosphere. This observable 
can be affected both by the proper increase of the p-p cross section at ultra- 
high energy (as in the case of the 1 / cosh profile) and also by the increase of 
the p-A cross section, due to the gluon saturation inside the target nucleus. 
The first possibility was discussed recently in |l7| . In our case, the p-p cross 
section fitted by the 1/cosh profile function is very close to the upper limit 



used in [l7j. However, as shown in Fig. 4, this mechanism is less effective 
than the GNSN scenario if the INGP description is applied. 

To calculate a realistic value _of (Xpav) we need a sophisticated simula- 



tion of the air- shower processes 
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20| involving all the exclusive cross 



sections. Here, just to get an idea on how the above increase of the cross 



section affects (X^ax) , we apply a simple toy model due to Heitler |2l| to 
estimate the deviation of (Xmax) from these calculations which are based on 
the Glauber description. Assuming such differences of (Xmax) can be identi- 
fied with the sum of differences of mean free paths calculated from the INGP 
and the GSNS scenario as 

AX^a.= Yl ^^(^^)' (53) 

where 

AX(E) = - (54) 

and {niair) is the average nuclear mass of the air, a^^ and 0"^, are cross 
sections calculated by the GNSN scenario and the INGP, respectively. The 
summation is done over the cascading steps while the energy of the leading 
particle is larger than the crossing energy Ecross from which the gluon sat- 
uration scenario overcomes the nucleon-nucleon Glauber cross picture. To 
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compare with the experimental data, we use the above toy model estimate 
for AXmax subtracting from the SIBYLL collaboration results for proton-air 



In Fig.5, we show the estimated (Xmax) values for three different pro- 
file functions (+ Gaussian, □ 1/cosh, x Woods-Saxson) together with the 
SIBYLL calculations for the proton-air and Fe-air simulations (dashed lines) 
and also the observed values extracted from the Pierre Auger Observatory 
experiment (black circles). 

5. Discussion and Perspectives 

In this paper, we have explored the idea of gluon saturation inside a 
nucleus in the high-energy limit and its effect on the proton-nucleus cross 
section. We show that if gluon saturation occurs in the surface region of the 
target nucleus, the proton-nucleus cross section starts to increase very rapidly 
as a function of the incident energy. Such a mechanism should eventually 
happen for an ultra-high energy scale, but the question is in which energy 
scale such a scenario starts to dominate. 

Such energy scale is determined by the form of the distribution of glu- 
ons near the surface area. If we assume that the small x gluon distribution 
in the nucleus follows that of the nucleon wave function inside the nucleus, 
the energy scale where the gluon saturation scenario starts to dominate the 
independent nucleon picture is around 10^^ — 10^^ eV. We note that dif- 
ferent profile functions give more or less the same energy scale once the 





(55) 
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proton-proton cross section is well fitted. It is very suggestive that the gluon 
saturation scenario inside the surface area seems consistent with the proton 
primary interpretation of the observed (Xmax) behavior in the Pierre Auger 
Laboratory experiment. 

As we see, the difference between the two pictures, INGP and GSNS, 
becomes very large at high energies. The reason for this is that, while in the 
INGP the effect of virtual gluons which bound the nucleon near the surface 
area is completely neglected, these gluons become dominant at high energies 
in the GSNS scenario. In a simple-minded argument, one might think that 
such an effect of nuclear binding must be negligible at high energies, since 
the ratio of the binding energy of a nucleon to the incident energy tends to 
zero. However, the situation may not be so simple. In the GSNS scenario, 
the density of virtual gluons, probably forming a kind of fractal fingers when 
penetrating into the vacuum among nucleons similar to the electric discharge 
pattern, becomes higher and higher at high energies, and eventually percolate 



everywhere even in 
condensate picture 



the nuc^ 
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ear surface region. According to the color glass 



25 



26|, such a scenario should happen at some 



energy scale, even at the lowest density region of the nuclear surface. 

Naturally, the energy scale depends on the precise form of the geometric 
gluon distribution inside the nucleus. If the distribution does not follow the 
probability distribution of nucleons but more sharp surface distribution, then 
the energy scale may shift to higher energy. In the example of Fig. 3 and Fig. 
4 the surface thickness parameter of the Woods-Saxon distribution is taken 
a little bit smaller than the usual value fitted to the nuclear distribution in 
nucleus {d ^ 0.6 fm) since this fit only applies for heavy nuclei {A > 40). 
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If we take d — 0.6 fm, the energy scale is lower by one order of magnitude. 
Therefore, the energy scale depends crucially on how the gluons starts to 
saturate in the nuclear surface region. Depending on this, the energy scale 
can be even lower than estimated here. To find a real energy scale where 
the gluon saturation occurs at the nuclear surface, further investigation on 
high-energy proton-nucleus or electron-nucleus collisions will be necessary. 

There exist ambiguities in the proton-proton cross section used in this 
paper (those of the model and extrapolation of the PDF, as well as those 
of experimental data) and these affect the precise value of the energy scale. 
However, the general conclusion of the present work does not change, as far 
as the energy dependence of the proton-proton cross section is fixed. 
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Tables 



Profile function 


Rd ifm) 


Rp ifm) 


C 


Qo (GeV) 


Gaussian 


0.602 


0.621 


2.8 


4.3 


1 / cosli 


0.61 


0.58 


3.0 


4.0 



Table 1: Parameters of the effective dipole description for the proton-proton 
cross section. 
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